Role of calcium and noise in the persistent activity of an isolated neuron 
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The activity of an isolated and auto-connected neuron is studied using Hodgkin-Huxley and 
Integrate-and-Fire frameworks. Main ingredients of the modehng are the auto-stimulating autap- 
tic current observed in experiments, with a spontaneous synaptic liberation noise and a calcium- 
dependent negative feedback mechanism. The distributions of inter-spikes intervals and burst du- 
rations are analytically calculated, and show a good agreement with experimental data. 



Understanding the mechanisms responsible for persis- 
tent activity and rhythm setthng is of central impor- 
tance in neuroscience. The persistence of the activity 
is the neuronal basis of the working memory [ij and 
brains rhythms ranging from about 0.1 to 200 Hz have 
been recorded in sleep, waking and pathological states 
[2|. Such behaviors are usually network properties; for 
instance oscillations are possibly built on the existence 
of different e.g. inhibitory and excitatory populations of 
neurons [3| ■ From a theoretical point of view, the large 
number of neurons in networks allows the use of self- 
consistent (mean-field) methods to determineproperties 
as the average spike emission frequency 0, IM Ifl S H ■ It 
was in particular found that a basis for persistent activity 
(with firing rate ~ 10 - 50 Hz) is the ability of NMDA 
synaptic channels to integrate afferent inputs with a slow 
decay time constant (~ 0.1 s) ^. 

Interestingly, persistence and rhythm settling have also 
been experimentally observed in systems made of an iso- 
lated and auto-connected excitatory neuron (autapse) [9| . 
Autapse are produced in vitro by grown excitatory neu- 
rons, extracted from rat embryos hippocampal, on cover- 
slips J9l linlllll |. Neurons normally develop up to 5 weeks 
and establish connections with themselves when no other 
neuron is nearby. The number of auto-connections in- 
creases with the age, as sketched in FigQ] Patch pipettes 
allows both electrical recording of the neuronal activity 
and current injection to trigger spikes. After a spike has 
been triggered, more than 2 week old autapses carry on 
spiking in a whole burst of activity. Records show that 
both the time interval between successive spikes (ISI) and 
the duration of the burst (BD) fluctuate. Surprisingly, 
while the number of auto-connections increase with the 
age, the average spike frequency decreases (around 20, 5, 
1 Hz for 2, 3, 4 week neurons respectively), see FigH 
The neuron therefore exhibits a negative rate-control 
feedback mechanism preventing runaway excitations due 
to the strong, and growing with the age, positive auto- 
stimulating current. All those experimental results are 
reported in [2j. 

This letter presents a theoretical study allowing a 
quantitative interpretation of the autapse activity. Per- 
sistence is due to the interplay of two currents evidenced 
in experiments !ll|: a small and slow postsynaptic com- 



FIG. 1: Sketch of the autaptic system (left panel) showing 
the increase in synaptic connections with age (2, 3 and 4 
weeks from top to down) and experimental (EX, middle) & 
theoretical (HH model, right) current time traces. Units are 
pA (EX) & /iA/cm^ (HH) vs. seconds. Experimental data 
show bursts of activity for 3 representative neurons among 20 
records. HH simulation show current dynamics for bursts of 
equivalent durations. Insets magnify two consecutive spikes 
chosen in the bursts, with typical ISI of 0.05, 0.2 and 1 s for 
2, 3 and 4 weeks respectively. The autaptic AMPA current 
(negative i.e. inward, immediately after the spike) is very 
small and not visible at 2 weeks, large and observable at 3 
weeks, huge at 4 weeks. 

ponent and a random spontaneous synaptic liberation. 
To model rhythm settling, a calcium dependent negative 
feedback is introduced. This feedback is at the origin 
of spike frequency adaptation under an external current, 
previously modeled Il2l | l3i and experimentally observed 
in the autaptic system |lll | . ISI and BD distributions are 
analytically calculated using an Integrate-and-Fire (IF) 
model, and compared to experimental data and numer- 
ical predictions from a detailed Hodgkin-Huxley (HH) 
model. 

Currents in the autaptic system are represented in 
FigUl Pre-synaptic (spike) currents are introduced in HH 
through standard Sodium and Potassium gating variables 
|l4.. 15j; their modeling in IF is discussed below. The ex- 
perimental characterization of the post-synaptic current 

has evidenced three auto-stimulating components that 
follow the spike (Fig.|21): 

1. A large amplitude excitatory AMPA component, 
I AM PA, entering just (5 ms) after the spike, and rapidly 



decaying. Its amplitude considerably increases with the 
age of the neuron, due to the growing number of con- 
nections (Fig. nj. The AMPA current is modeled, in 
HH, through an effective exponentially decreasing con- 
ductance [15|. We stress that the AMPA current arrival 
falls within the neuron refractory period and thus can- 
not by itself trigger a new spike; it however leads to a 
membrane potential depolarization (up to mV), giving 
a bump in the flank of the spike in voltage records [2|. 
The AMPA depolarization has two consequences. It first 
slows down Na and K gating variable resets, increasing 
the refractory period. Secondly, it allows more calcium 
to enter the cell via voltage-gated channels, as soon as 
membrane potential exceeds -20 mV. In IF, AMPA, Na 
and K currents are altogether accounted for by the cal- 
cium Cw^^ entering at each spike, and the refractory pe- 
riod tr- These two parameters (which depend on the age 
of the autapse) and the threshold 9 for spike fi ring are 
determined from the numerical analysis of HH [l^l, see 
Fig. El and [13. 

2. A small amplitude and slow decaying component, 
Id: due to NMDA and ICAN conductances 0. Id de- 
polarizes the membrane potential V during a burst to 
~ —55 mV; after the burst halts, V decays to the rest 
value Vi = —60 mV with a time constant to ranging from 
300 ms in young cells (2 weeks) to 1 s in mature cells (3 
and 4 weeks). In both HH and IF, lo is modeled as an ex- 
ponentially decaying current, loit) = I^ e^'^*^*'^'^*^/*°, 
released with a delay St = 5 ms after the spike occurring 
at time i^; the amplitude is set to I^) = 2 mA/cm^, and 
corresponds to a depolarization of the membrane poten- 
tial of 5 mV. 

3. A spontaneous and random synaptic release, called 
miniatures. Is = ggV s{t). Both amplitude and fre- 
quency of the miniatures s{t) are stochastic |2|. The 
time interval between miniatures during a burst has been 
fixed to its average value 6s ~ 20 and 10 ms for 2-3 and 
4 week neurons respectively. 6s increases to its (much 
larger) rest value ~ 0.1 s in about 5 s after the end of the 
burst. Release times are hence discrete and multiple of 
6s, which makes the IF model mathematically tractable 
(taking into account the stochasticity in times between 
miniatures does not significantly affect the outcome |l6|). 
The distribution P of the released amplitude a is Poisso- 
nian; the mean m and the conductance gs are fitted from 
experiments jl3| , and given in caption of Fig. El After re- 
lease at time i x 6s, the miniature exponentially decays, 
s(t) = {s- +a{i)) e-(*-"'-)/*= with t^ = 5 ms (< 6s) 
from experiments. Parameter s~ — 77i/(e^'*'' — 1) is the 
average residual amplitude before the release. 

In addition, an inhibitory feedback due to the presence 
of a Calcium-dependent After Hyperpolarizing current, 
Iahp — gAHP j^ {^ ~ Vk), is introduced. Such a com- 
ponent is modeled as in [l^l for HH and |j,ll3 for IF. Pa- 
rameters values are gAHP = 5 mS/cm^, Vk = — 80 mV, 
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FIG. 2: Currents in the auptatic neuron model, with line 
widths proportional to amplitudes. Presynaptic currents Inu 
and Ik accompany a spike. Postsynaptic components are, 
see text: I am pa (thick line), responsible for membrane de- 
polarization and more Ca entering the cell (wiggly line); Id, 
a slow depolarizing current (thin line); Is, due to enhanced 
random spontaneous miniatures following Ca release (dashed 
line). The negative feedback is accounted for by a potas- 
sium hyperpolarizing current, Iahp, with Ca-dependent am- 
plitude. 
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FIG. 3: Phase diagram for the 
{to, Tags) plane for gi — 0.1,9 



persistence of activity in 
= -53.5 mV (left) and 



gi = 0.2, 9 — —52 mV (right). Bars symbolize ranges of exper- 
imentally observed values Q • Parameters are: Ca''''=0.043 (2 
week, full), 0.115 (3 week, dotted) and 0.55 fiM (4 week neu- 
ron, dashed line). Activity persists above the lines, is absent 
below. Left: parameters, indicated by the star in (0.3,0.004) 
with gs = 0.002, m = 2, for a 2 week cell with stationary 
ISI of 0.07 s. Right: stars lies slightly below threshold at 
(1,0.02) ~gs = 0.0032, m = 6.5 (3 weeks)- and (1,0.022) - 
gs = 0.0032, m — 7 {4 weeks). Positive fluctuations of the 
noise are needed for the activity to persist. Inset: magnifica- 
tion of the star-surrounding region; arrows indicate a positive 
standard deviation reaching (3 weeks: narrow line on the left) 
or (4 weeks: bold line on right) crossing the critical line. Av- 
erage ISI are 0.2, 0.75 s for 3, 4 weeks respectively. 

kd = 30 ^M. 

In the IF model, the dynamics of membrane potential 
V and calcium concentration Ca obey 



C^^-Ii-lAHp{Ca) 
at 



Id - Is 



^^Ca^py^Sit-t,,)^^ 

at ^-^ tca 



(1) 

(2) 



where Ii — gi (V — Vi) is the leak current, C — 1 mF/cm^ 
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FIG. 4; From top to down, 2, 3 & 4 week neurons. Left; 
Calcium levels vs. time during a burst for random -HH, nu- 
merical, full line- and fixed amplitude -IF, theory, circles- 
miniatures. Jumps coincide with spikes. Right: IF -theory, 
full- and HH -numerical, dotted line- distributions of calcium 
right after spike emission. Distribution is very peaked for the 
2 week neuron around Ca — 0.26 (IF) and 0.28 (HH), and 
much wider due to noise-driven activity in older cells. 

the capacity, and t„ the emission time of the n*'* spike. A 
spike is fired when V > 0, and the refractory period tr is 
over. Large amphtude synaptic miniatures {Is) that add 
to the slow component {Id) can trigger a spike, provided 
the calcium-based inhibition {Iahp) is sufficiently weak. 
This mechanism is made possible by the fact that the 
slow component and high-frequency miniatures persist 
over a time period larger than (or, at least, equal to) the 
calcium decay constant tca = 0.33 s. 

This few hundred ms time-scale is well separated from 
the integration time, r,„ — C /gi = 5 — 10 ms for 
gi = 0.1 — 0.2, and miniature decay time, t^ — h ms . 
Consequently, once time t after the last spike emission is 
expressed in terms of miniature intervals, t — i x 5s + t, 
where i{> 0) is integer- valued and r continuously ranges 
in [0; 5s\, Iahp and Id do not vary with r 18]. From |(2Jl, 
the calcium level Ca is a function of its value Ca^ right 
after the previous spike and the interval i: Ca[Ca-, i] = 
Ca- e^''-^=/*-ca _ -y^g then integrate (Q) analytically over 
T at fixed lD{i),lAHp{i),<^{i) and Ca-. The resulting 
expression for the maximum of the membrane potential 
over r, Va/(z, <y{i), Ca-) [ig, is then compared with 6. i 
is increased by one if no spike is fired; otherwise, it is re- 
set to and the calcium level is increased by Ca^^. Spike 
emissions (Fig.QJ are therefore reflected by jumps in the 
Ca level, as illustrated in Fig. ^ 

Insights about the values of parameters allowing ac- 
tivity to be persistent can be obtained when consid- 
ering first non-stochastic miniatures with fixed ampli- 
tude m. In this case, a stationary dynamics correspond- 
ing to a regular burst with infinite duration can set- 
tle down. The values for the calcium after each spike, 
Ca*, and ISI, i* 5s, are such that the calcium decrease 



between two spikes balances the increase after a spike, 
Ca* = Ca^P/{l - e-^**»/*ca). Inserting this result into 
the IF condition Va/ («*, ^ti, Ca*) = 9 derived above yields 
a self-consistent equation for the stationary ISI value i* . 
The existence or absence of solution to this equation in- 
dicates whether a persistent activity can be sustained or 
not [l9j- The phase diagram is shown in Fig.O for two 
values of the neuron conductance, gi = 0.1,0.2, in the 
experimentally measured range. 

The presence of noise in miniatures remarkably en- 
riches this picture. Experimental features of ISI (aver- 
age value in Fig. ^ small fluctuations for 2 week neurons 
and wide distribution for 3-4 weeks in Fig. O are well 
reproduced by parameters shown in Fig. |3| The system 
representative parameters (star) lies above the critical 
persistence activity line for 2 week cells, and slightly be- 
low for 3-4 weeks. In the latter case, miniatures show- 
ing positive fluctuations occasionally make representative 
parameters cross the phase border, and the neuron gener- 
ates spike (noise-driven activity). On the contrary, in the 
former case, activity would persist forever in the absence 
of noise. 

Finite burst duration results from a two-fold mecha- 
nism. On the one hand, if miniature amplitudes a{i) 
happen to be lower than the minimal value that can trig- 
ger a spike in the vicinity of i* , no spike will be fired. 
Later on, as the slow component decreases, spike firing 
becomes less and less likely. On the other hand, as a 
result of positive fluctuations in the its, a very short ISI 
may arise, making calcium increase and further spike fir- 
ing unlikely. 

The above results are illustrated in Fig. ^ for both 
HH/IF models, and noiseless/noisy dynamics. Note that, 
for 2 weeks neurons, persistent activity is possible in the 
absence of miniatures too (m = 0) through the slow com- 
ponent only for 2 week cells. A burst, in the noisy dy- 
namics, halts after a positive increase in calcium due to 
contiguous spikes with short ISI. For 3 and 4 week cells 
the noiseless IF dynamics does not reach stationarity, 
while the noisy dynamics results in finite and random 
duration bursts. 

To calculate the stationary distribution of calcium, 
ISI, BD in presence of noise, we consider the prob- 
ability T{i\Ca-) to fire with an ISI equal to i5s 
given the calcium value Ca- after the previous spike: 
mCa-) = T{7=\pU\Ca-) (1 -p(z|Ca_)); p(j|Ca_) = 

Sfc^o " P{^) is the probability not to fire at in- 

terval j where Sm{j\Ca-) is the minimal amplitude of 
noise able to trigger a spike at interval j, P is the 
Poissonian with average m. The transition matrix be- 
tween two successive calcium values read T{Ca\Ca-) = 
Y.^f{i\Ca-)5{Ca - Ca[Ca-,i] - Ca^P). The station- 
ary calcium distribution, Q{Ca), is the maximal eigen- 
vector of T. We have discretized the calcium interval, 
and used Kellogs iterative projection method to diag- 



onalize T with the result shown in Fig. 0] The ISI 
distribution can be obtained from the matrix product 
of T{i\Ca) and Q, and is shown in Fig. [S] It is very 
peaked around 0.05 (EX) - 0.06 (IF) s for 2 week cells, 
and spread out for 3 and 4 week neurons, with median 
0.2, 0.75 s. respectively, in good agreement with exper- 
iments. To calculate the BD distribution, we define the 
generating function for the probability of an interval i Ss 
between spikes, G{x) — J2i Q{i)x^ ■ The coefficient of x'^ 
in G{x)^, denoted by [x'']G{xy^ is the probabihty that 
a burst with n spikes has duration k5s- Summing over 
n, we get the probability that a burst has duration k5s^ 
[x''] G{x)/{\ - G{x)). The resulting BD distribution is 
shown in Fig. El it decays as e~^l^° where ko is the root of 
(^(•gi/feo~) — \ ^}jg agreement with experiments is good, 
even if comparison suffers from limited data (9, 13, 11 
bursts for 2, 3, 4 weeks); it is excellent with HH simu- 
lation, which makes the analytical study of IF quite at- 
tractive despite the approximations done (discretization 
of time) . 

In conclusion, this letter proposes a possible mecha- 
nism that accounts for the persistence of activity in an 
isolated autapse. The slow current In sets the neuron 
in a depolarized 'up' state, also observed in various os- 
cillatory networks |3 where it results from e.g. afferent 
inputs from other neurons. For 2 week neurons this up 
state can, by itself, sustain activity while, for 3 and 4 
week, activity requires noisy synaptic liberation. Noise 
is at the same time responsible for persistent dynamics, 
robust with respect to changes of parameters [l9j, and 
the finite duration of bursts. The noise-driven mecha- 
nism for burst halting presented here differs from mech- 
anisms based on a slow activity dependent depression 
P (either due to a modulation of cellular excitability, 
or to synaptic depletion [23), and is supported by the 
observed large variability in burst durations. In addi- 
tion, calcium-dependent spike frequency adaptation can 
explain the observed pattern of activity as the increase 
of ISI with the age. Finally, it would be interesting to 
extend the present study to the case of a network com- 
posed of a small number of similar neurons for which 
experimental data are available. 
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